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Abstract. Gamma ray bursts (GRBs) belong to the most energetic 
events in the Universe. Recently, the extragalactic nature of these sources 
has been confirmed with the discovery of several host galaxies (HGs) and 
the measurement of their redshifts. To explain the origin of GRBs various 
models have been proposed, among which the coalescence of compact 
objects and the collapsar scenarios are the most representative, being the 
collapsar model one of the most accepted to explain the long duration 
GRBs. A natural consequence of this model is that the GRBs would 
trace the star formation rate (SFR) of their HGs. In this contributed 
paper we present preliminary results of the development of a Montecarlo- 
based code for collapsar event formation which is coupled to chemical- 
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cosmological simulations aiming at studying the properties of HGs in a 
hierarchical scenario. 

Resumen. Las explosiones de radiacion gamma (ERG) se encuentran 

entre los eventos mas energeticos del Universe. Recientemente, la natu- 
raleza extragalactica de estas fuentes fue confirmada con el descubrim- 
iento de varias de sus galaxias huesped y la medicion de sus corrimientos 
al rojo. Para cxplicar su origen, varies modelos han sido propuestos, entre 
los cuales las colisiones de objetos compactos en sistemas binaries y los 
llamados collapsars son los mas representatives, siendo este ultimo med- 
ele uno de los mas aceptados para explicar las ERG dc larga duracion. 
Una consecuencia natural del mismo es que las ERG trazan'an la tasa de 
formacion estelar de sus respectivas galaxias. En esta comunicacion se 
presentan rcsultados prcliminares del desarrollo dc un codigo Montccarlo 
para la formacion de collapsars en simulaciones qufmico-cosmologicas, con 
el fin de estudiar las propiedades de sus galaxias huesped en un escenario 
jerarquico. 

1. Introduction 

The gamma ray bursts (GRBs) are the most energetic electromagnetic events 
in the Universe (e.g. Piran et al. 2000). One of the preferred models to ex- 
plain the long duration (i.e., > 2 s) GRBs is the star core collapse into a black 
hole produced in a supernova type SNIb/c explosion (Mac Fayden, Woosley & 
Heger 2001). This mechanism of GRB formation is also known as the coUapsar 
scenario. This model is linked to the evolution of massive stars which normally 
have a mean lifetime of several million years implying that the typical lifetime of 
a GRB progenitor system is negligible in cosmological terms. So, it turns to be 
natural to consider these events as possible tracers of the cosmic star formation 
history up to high redshifts. In particular, the GRBs would permit to obtain 
information about star formation regions in galactic systems with different levels 
of evolution. 

Recently, Courty et al. (2004) made use of structure formation simulations 
in order to identify galactic populations in the simulated sample capable to 
reproduce the observational features of the observed HGs. These galaxies show 
a trend to be bluer and sub-luminous (e.g. Le Floc'h et al. 2003). Their analysis 
seems to confirm the connection with the star formation rate (SFR) if the GRBs 
events are formed in galaxies with high star formation efficiency. In this work we 
develop a GRB event generator based on the collapsar model for the progenitors 
and study the properties of the HGs in cosmological simulations. 



2. GRB Montecarlo code and discussion 

The GRB event algorithm has been designed to work coupled to numerical 
simulations in a cosmological framework. We used the simulated galaxy catalog 
constructed by De Rossi et al. (2005) from simulations performed with the 
chemo-dynamical code of Scannapieco et al. (2005) developed within GADGET-2 
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Figure 1. Left: Cosmic SFR density for the simulated box of 10 Mpc 
comoving size as a function of redshift. Right: Simulated proba- 
bility for GRBs occurrence as a function of redshift. 



(Springel & Hernquist 2002). The gas component is followed using the Smoothed 
Particle Hydrodynamics technique while the dark matter component is followed 
using a Trcc-Particlc Mesh method. This chcmo-dynamical code describes the 
enrichment of the interstellar medium by SNII and SNIa supernova explosions. 
The cosmological model assumed is defined by the following set of cosmological 
parameters: Qm = 0.3, Qa = 0.7, Jib = 0.04 and Hq = 100 h km s~^ Mpc~^ with 
h = 0.7. The simulations represent a typical region of 10 Mpc comoving size 
with 2 X 80^ particles which translates into an initial mass resolution of 2 x 10^ Mq 
and 2 x 10^ Mq h^^ for the gas and dark matter particles, respectively. The 
catalog of simulated galaxies provides information on the gas, stellar and dark 
matter components for the redshift range 2; = [0, 3]. A Montecarlo-based code 
was developed to generate GRBs events in each simulated galaxy of the catalog. 
We assumed Poissonian statistics to emulate the probability distribution of the 
GRBs events. The GRB generator selects young stars as possible candidates, 
where SNIb/c events may occur, and consequently, collapsar events can develope. 
For that, we adopted a cut-off in stellar ages of tc = 10^ yr. A median rate of 
collapsar events (i?coii) consistent with that estimated by Fryer et al. (1999), 
-Rcoii ~ 10 — 1000 Myr~^ Galaxy"^ where Galaxy represents a typical galaxy, 
was assumed. In particular, in this work we used i?coll = 100 Myr~^ Galaxy"^ 
and a typical galaxy mass of 10^^ Mq. A total of 500 Montecarlo realizations 
was perfomed for every selected stellar particle in each simulated galaxy in the 
catalog. Our preliminary results can be seen in Figures 1 and 2. 
Figure 1 shows the simulated comoving SFR density (psFR) left panel) and the 
simulated probability of GRB event occurrence (right panel) as a function of 
redshift. As it can be seen from Figure 1, assuming the collapsar scenario, 
where GRBs originate from massive stars, produces the expected behaviour of 
GRBs being good tracers of the cosmic SFR history. Figure 2 shows the SFR 
efficiency (csfr, defined as the SFR normalized to the total stellar mass at each 
analysed redshift) versus circular velocities (Kjpt) of the HGs from z = 3. It can 
be seen from this figure that the SFR efficiency anticorrelates with V^pt , so that 
the higher the efficiency, the smaller the systems (i.e. slower rotating systems). 
Note that, at all analyzed z, more than 50 per cent of the simulated HGs have 
V^pt < 100 km s~^. According to the Tully-Fisher relation, this result implies 
that high SFR efficiencies are associated to sub-luminous systems in agreement 
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Figure 2. SFR efficiency versus circular velocity for the simulated 
HGs sample from z = 3 to z = 0.1. 

with observations. Note also that the SFR efficiency increases with increasing 
redshifts and that for z > 1, all HGs seem to have similar values. 
Interesting improvements for the future would be to include an analysis based 
on the simulated colours and metallicities of the simulated HGs. Numerical 
simulations with higher resolution are also being analysed in order to confirm 
these trends. 
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